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SPREAD SPECTRUM MULTIPATH COMBINATION 



CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims priority from provisional application Serial No. 
60/152,974, filed 09/09/99. The following patent applications disclose related 

subject matter: Serial Nos. 08/ .., filed .... (Tl — ). These applications have a 

common assignee with the present application. 

BACKGROUND OF THE INVENTION 

The present invention relates to digital communications, and more 
particularly to multipath received signal combination systems and methods. 

Spread spectrum wireless communications utilize a radio frequency 
bandwidth greater than the minimum bandwidth required for the transmitted data 
rate, but many users may simultaneously occupy the bandwidth. Each of the 
users has a pseudo-random code for "spreading" information to encode it and for 
"despreading" (by correlation) received spread spectrum signals and recovery of 
information. Figure 2 shows a system block diagram, and Figure 3 illustrates a 
QPSK (quadrature phase-shift keying) encoder with pseudo-noise spreading. 
Such multiple access typically appears under the name of code division multiple 
access (CDMA). The pseudo-random code may be an orthogonal (Walsh) code, 
a pseudo-noise (PN) code, a Gold code, or combinations (modulo-2 additions) of 
such codes. After despreading the received signal at the correct time instant, 
the user recovers the corresponding information while other users' interfering 
signals appear noise-like. For example, the interim standard IS-95 for such 
CDMA communications employs channels of 1 .25 MHz bandwidth and a 
pseudo-random code pulse (chip) interval T c of 0.8138 microsecond with a 
transmitted symbol (bit) lasting 64 chips. The recent wideband CDMA (WCDMA) 
proposal employs a 3.84 MHz bandwidth and the CDMA code length applied to 
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each information symbol may vary from 4 chips to 256 chips. Indeed, UMTS 
(universal mobile telecommunications system) approach UTRA (UMTS terrestrial 
radio access) provides a spread spectrum cellular air interface with both FDD 
(frequency division duplex) and TDD (time division duplex) modes of operation. 
UTRA currently employs 10 ms duration frames partitioned into 15 time slots with 
each time slot consisting of 2560 chips (T c = 0.26 microsecond). 

The CDMA code for each user is typically produced as the modulo-2 
addition of a Walsh code with a pseudo-random code (two pseudo-random 
codes for QPSK modulation) to improve the noise-like nature of the resulting 
signal. A cellular system as illustrated in Figure 4 could employ IS-95 or 
WCDMA for the air interface between the base station and the mobile user 
station. 

A spread spectrum receiver synchronizes with the transmitter by code 
acquisition followed by code tracking. Code acquisition performs an initial 
search to bring the phase of the receiver's local code generator to within typically 
a half chip of the transmitter's, and code tracking maintains fine alignment of chip 
boundaries of the incoming and locally generated codes. Conventional code 
tracking utilizes a delay-lock loop (DLL) or a tau-dither loop (TDL), both of which 
are based on the well-known early-late gate principle. 

The air interface leads to multipath reception, so a RAKE receiver has 
individual demodulators (fingers) tracking separate paths and combines the 
finger results to improve signal-to-noise ratio (SNR). The combining may use a 
method such as the maximal ratio combining (MRC) in which the individual 
detected signals in the fingers are synchronized and weighted according to their 
signal strengths or SNRs and summed to provide the decoding statistic as 
illustrated in Figures 5a-5b. That is, a RAKE receiver typically has a number of 
DLL or TDL code tracking loops together with control circuitry for assigning 
tracking units to the strongest received paths. Also, an antenna array could be 
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used for directionality by phasing the combined signals from the antennas. 

However, known multipath combination methods have problems including 
robustness. 

SUMMARY OF THE INVENTION 

The present invention provides a multipath received signal combination 
covariance matrix eigenvector components used for relative path weightings. 
Linear approximations of channel response provide models for covariance matrix 
analysis. 

This has the advantages of a robust maximal ratio combination. Also the 
technique can be blind; that is, it does not require pilot symbols. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The drawings are heuristic for clarity. 

Figures 1a-1c show a preferred embodiment detector, linear 
approximation, and flow diagram for multipath combining. 

Figure 2 shows a CDMA communications in block diagram. 

Figure 3 illustrates a.CDMA encoder. 

Figure 4 indicates cellular communications. 

Figures 5a-5b illustrate RAKE receiver fingers and antenna array. 



TI-29648 Page 3 



DESCRIPTION OF THE PREFERRED EMBODIMENTS 



Overview 

Preferred embodiment spread spectrum communications systems use 
preferred embodiment blind (without synchronization) maximal ratio combining 
(MRC) methods with multipath detectors. Preferred embodiment MRC methods 
estimate time varying multipath channels using eigenvectors associated with 
maximal eigenvalues of matrices correlating individual path outputs. Local linear 
approximations of channel parameter variation overtime simplify the matrix 
analysis. Figure 1a illustrates a preferred embodiment RAKE detector which 
may apply the preferred embodiment MRC methods. Figure 1b suggests the 
channel parameter linear approximation, and Figure 1c is a flow diagram. In 
contrast, the known MRC methods estimate each path separately and then 
combine with weightings proportional to path output power. 

The preferred embodiment MRC methods include updating channel 
estimates both on a symbol-by-symbol basis and on a slot (sequence of 
symbols) basis. For an updating on a symbol-by-symbol basis, the path 
covariance analysis averages over output sequences from each path where the 
output sequences have length 4K+1 and are centered about the symbol interval 
for updating; thus each updating requires new covariance matrices. Typically, 
for WCDMA take K in the range of roughly 20 to 40, depending upon the Doppler 
rate and the symbol rate (32 KSPS channel). For an updating on a slot basis, 
the channel estimates for all of the 4K+1 symbol intervals in the slot derive from 
the a single set of covariance matrices over the length 4K+1 output sequences. 
This amounts to a linear approximation of the channel over the slot. 

Interim standard IS-95 and wideband CDMA (WCDMA) proposals include 
a pilot signal which is much stronger than the data signals and can be used for 
preferred embodiment overall channel estimation. 
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In preferred embodiment communications systems base stations and 
mobile users could each include one or more digital signal processors (DSP's) 
and/or other programmable devices with stored programs for performance of the 
signal processing of the preferred embodiment synchronization methods. 
Alternatively, specialized circuitry could be used. The base stations and mobile 
users may also contain analog integrated circuits for amplification of inputs to or 
outputs from antennas and conversion between analog and digital; and these 
analog and processor circuits may be integrated on a single die.. The stored 
programs may, for example, be in ROM or flash EEPROM integrated with the 
processor or external. The antennas may be parts of receivers with multiple 
finger RAKE detectors for each user's signals. Exemplary DSP cores could be in 
the TMS320C6x or TMS320C5x families from Texas Instruments. 

Symbol-by-symbol update preferred embodiments 

Preferred embodiments estimate the multiple paths of a received signal 
and perform MRC for detection of a transmitted symbol s(0) from a sequence of 
transmitted symbols s(k). The received baseband signal at a RAKE detector can 
be expressed as the real part of 

r(t) = E k I m h m (k) s(k) c(t-kT-T m (k)) + noise(t) 
where T denotes the time interval of a symbol, c() denotes the pseudo-random 
code of values ±1 (and equals 0 outside of the interval [0,T]), and m indexes the 
multipath with relative delay x m (k) and channel gain and phase h m (k) during 
symbol s(k) propagation. Presume normalized symbols: |s(k)| 2 = 1. Each RAKE 
detector finger tracks and despreads (correlates) one of the multipaths using cQ 
to despread and synchronize. Let y m (k) = h m (k) s(k) + v m (k) denote the mth 
finger output where v m (k) models noise along the mth path. Indeed, take the 
v m (k) to be independent, normally distributed, and with standard deviation c. 
Then MRC weights the outputs y m (k) according to their power and phase, so 
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ideally the MRC weightings are proportional to h m (k)* and the MRC decoding 
statistic with 0 the overall phase would be 

e ie I m y m (k)h m (k)* 
The expectation of y m (k) over a sequence of ks is approximately 0 
because s(k) is randomly +1, +j (assuming QPSK modulation) and varies more 
rapidly than h m (k); thus consider the MxM covariance matrix C of the 
synchronized finger outputs with components given by 

C m . n = (I k y m (k)y n (k)*)/(4K+1) 
where the k sum is over the range -2K < k < 2K with K on the order of 20 to 40. 
Expressed in matrix form with y(k) the M-dimensional vector of synchronized 
outputs and (.) H denoting Hermitian conjugate (adjoint): 

C = I k y(k)y H (k)/(4K+1) 
Locally, h m (k) can be reasonably approximated by low order polynomials in k. 
Thus consider constant and linear approximations. 

1 . First consider a time-invariant channel; that is, take h m (k) = h m for all k in the 
range -2K to 2K. With this presumption and using the fact that s(k) is randomly 
±1 , ±j, the covariance matrix then becomes: 

C = hh H + <j 2 I 

where h is the M-component vector with mth component equal to h m . This 
essentially is a spectral resolution for C and shows that C has one eigenvalue 
equal to |h| 2 + a 2 which is associated with the normalized eigenvector h/|h| plus 
an (M-1 )-order eigenvalue equal to a 2 which is associated with M-1 (generalized) 
eigenvectors orthogonal to h. Thus the preferred embodiment proceeds as 
follows: 

(1) form covariance matrix C as averages over sequences of finger 
outputs. 

(2) find A, MAX , the largest eigenvalue of C, and its associated normalized 
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eigenvector, w, and 

(3) define the MRC (up to an overall phase) decoding statistic to be 

2m w* m y m (0) 

(4) determine the overall phase of this MRC statistic by any channel 
estimation method, such as taking fourth powers to eliminate the symbol phase 
(presumed QPSK modulation) and using a phase-locked loop. 

(5) decode the symbol s(0) from the resultant phase-adjusted MRC 
statistic. 

2. Next, consider the case of h(k) varying linearly with time. In particular, 
presume that for k in the range -2K to 2K: 

h(k) = a + kp 

where a and p are (unknown) constant M-dimensional vectors. Figure 1b 
heuristically illustrates the approximation for a component h m (k) with the broken 
horizontal line representing ot m and the slope of the diagonal line p m . So 
approximate h(k) and thereby determine the MRC weightings from estimates for 
a and p. Thus first define covariance matrices C 0 and as the sums over the k 
intervals -K < k < -1 plus 1 < k < K and -2K < k < -(K+1 ) plus K+1 < k < 2K, 
respectively: 

C 0 = ^ k=liK y(k)y H (k)/2K 

C 1 = I ±k=K+ i l2 Ky(k)y H (k)/2K 
Figure 1b illustrates the two portions of the -2Kto 2K symbol interval slot used to 
define the two covariance matrices. Then invoking the presumed linear form of 
h(k) yields: 

C 0 = aa H + f 0 (K) PP H + r 0 + % 

C 1 = aa H + f^K) pp H + r, + % 
where f 0 (K) = (2K 2 +3K+1)/6 and f,(K) = (14K 2 +9K+1)/6 are the corresponding 
sums of k 2 over the two k ranges, each divided by 2K, and the noise terms r 0 , 
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*F 0 , r 1f and ^ arise from the y(k)-noise cross terms and noise squared terms in 
the two sums of y(k)y H (k) over k. These variables have the following properties: 
mean(r 0 ) = meanfl^) = 0; mean(T 0 ) = meanOFJ = 0; var(r 0 ) = varfr,) = 
a 2 ||a|| 2 l/2K; and varOF 0 ) = var^FJ = cr 2 l/2K; thus consistent estimates can be 
obtained. 

Estimate a and p essentially by ignoring the noise terms in the foregoing 
pair of linear equations and solving for ococ H and pp H followed by eigenvector 
analysis. In particular, first define the estimates 

A= (f^KJCo-fo^C^KJ-foCK)) 
B= (-Co + C^MKJ-WK)) 
Thus A estimates aa H + noise terms and B estimates pp H + noise terms. As 
previously described, the largest eigenvalue of A equals ||a|| 2 + noise variance 
and has a corresponding normalized eigenvector which estimates a/||a|[ up to an 
overall phase. Similarly, the largest eigenvalue of B equals ||p|| 2 + noise 
variance and has a corresponding normalized eigenvector which estimates p/||p|| 
up to an overall phase. 

In short, the preferred embodiment method of MRC for symbol s(0) uses 
the linearized model of h(k) and includes the following steps: 

(1 ) compute the covariance matrices C 0 and from the vector y(k) of 
synchronized detected RAKE finger signals for symbol index k in the range -2K 
to 2K; k = 0 corresponds to the symbol of interest for the MRC. That is, 

Co = S± k= i, K y(k)y H (k)/2K 
C 1 = S± k =K + i,2Ky(k)y H (k)/2K 

(2) generate the estimate A from the covariance matrices C 0 and by 

A= (WJCo-fofKJC^KJ-WK)) 
where f 0 (K) = (2K 2 +3K+1)/6 and f^K) = (14K 2 +9K+1)/6. 

(3) find the largest eigenvalue, X A , and corresponding normalized 
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eigenvector, w A , of the estimate A. 

(4) define the MRC (up to an overall phase) decoding statistic to be 

S m w\ m y m (0) 

(5) determine the overall phase of this MRC statistic by any channel 
estimation method, and 

(6) decode the symbol s(0) from the resultant phase-adjusted MRC 
statistic. 

Update per slot preferred embodiments 

The preferred embodiment method of MRC for all ks in a slot of k values 
from -2K to +2K (such as illustrated in Figure 1b) again uses the linearized model 
of h(k) and includes the following steps: 

(1 ) compute the covariance matrices C 0 and C, from the vector y(k) of 
synchronized detected RAKE finger signals for symbol index k in the range -2K 
to+2K. That is, 

C 0 = E± k=1iK y(k)y H (k)/2K 

c 1 = WyWy"(P 

(2) generate the estimates A and B from the covariance matrices C 0 and 

Ci by 

A= (MKJCo-foOQC^^-WK)) 

where f 0 (K) = (2K 2 +3K+1)/6 and f^K) = (14K 2 +9K+1)/6. 

(3) find the largest eigenvalues, X A and X Bi and corresponding normalized 
eigenvectors, w A and w B , of A and B, respectively. The eigenvalues are real due 
to all of the matrices y(k)y H (k) being Hermitian. 

(4) estimate the channel for k in the slot as h'(k) = Vx, A w A + kA B w B . 
Because w A and w B are only determined up to a phase, use the k=0 case as 
previously described to determine the overall phase for w A . When s(k) is QPSK, 
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we can compute I ±k=1 , K (y T (k)y(k)) 2 /2K and Z+^i^ (y T (k)y(k)) 2 /2K. The 
expressions obtained will be functions of w A T w A , w B T w Bl X A , X B and the residual 
phases on w A and w B . The residual phases can be determined from the resulting 
equations since we know all the other quantities. In case s(k) is BPSK we only 
need to use compute Z± k=1 , K y T (k)y(k)/2K and E± k=K+ i.2K y T (k)y(k)/2K. A 
combination of this and the phase estimation method for k=0 can also be used to 
estimate the residual phases. Alternatively, pilots from the transmitter could be 
used for phase determinations. 

(5) use the components of h'(k) for relative weightings in the MRC of y(k) 
for each k in the slot (k=0,K was already done in step (4)). That is, take the MRC 
statistic for each k in the slot to be the sum 

S m y m (k)h' m (k)* 

(6) for each k decode the symbol s(k) from the resultant phase-adjusted 
MRC statistics. 

Alternative preferred embodiments 

Alternative preferred embodiment quadratic (or higher order) 
approximations for the channel presume h m (k) = <x m + kp m +k 2 y m and thus divide 
the sequence of received symbols into three sections and form three covariance 
matrices to solve for estimates of aoc H , pp H + ay H + ya H , and yy H followed by 
eigenvalue and eigenvector analysis to estimate a, p, and y and thus h(k) for 
use in the MRC. Of course, taking k in the interval -3K to 3K makes the notation 
simpler and there are three coefficient factors f 0 (K), f^K), and f 2 (K) as second 
order polynomials and three more coefficient factors as fourth order polynomials. 

Modifications 

The preferred embodiments can be modified in various ways while 
retaining the features of a finding relatively weightings for multipath combining by 
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eigenvector components of path covariance matrices. 

For example, .instead of updating the MRC coefficients every symbol one 
can update them every L'th symbol and interpolate or use the same estimates 
for the symbols in between. 

Similarly instead of every slot one can compute every half slot or two slots 
and interpolate for the values in between. 

Other techniques for phase estimation for e.g., decision directed phase 
estimation. 

Other weightings in the multipath combinations could be used; this could 
involve (fractional) powers of the channel estimation. 
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CLAIMS (29648) 



What is claimed is: 

1 . A method of multipath combining, comprising: 

(a) forming at least one matrix of covariances of multipath inputs; 

(b) finding an eigenvector of said matrix; and 

(c) combining said multipath inputs relatively weighted according to the 
components of said eigenvector. 

2. The method of claim 1, wherein: 

(a) said eigenvector is associated with a maximal eigenvalue of said 

matrix. 

3. The method of claim 1 , wherein: 

(a) said weightings have magnitudes proportional to the squared 
magnitudes of said eigenvector components and have phases proportional to the 
phases of said eigenvector phases. 

4. A method of multipath combining, comprising: 

(a) forming a first matrix of covariances of multipath inputs over a first 
range and a second matrix of covariances of multipath inputs over a second 
range; 

(b) forming an estimation matrix from said first and second matrices; 

(c) finding an eigenvector of said estimation matrix; and 

(c) combining said multipath inputs relatively weighted according to the 
components of said eigenvector. 
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5. The method of claim 4, wherein: 

(a) said eigenvector is associated with a maximal eigenvalue of said 
estimation matrix. 

6. The method of claim 4, wherein: 

(a) said weightings have magnitudes proportional to the squared 
magnitudes of said eigenvector components and have phases proportional to the 
phases of said eigenvector phases. 

7. The method of claim 4, further comprising: 

(a) forming a second estimation matrix from said first and second 
matrices; 

(b) finding a second eigenvector of said second estimation matrix; and 

(c) wherein said combining said multipath inputs relatively weighted 
according to the components of said eigenvector includes relatively weighted 
also according to the components of said second eigenvector. 

8. A multipath receiver, comprising: 

(a) a plurality of detectors, each detector detecting a path of a multipath 

input; 

(b) first circuitry coupled to said detectors and connected to form a first 
matrix of covariances of outputs of said detectors over a first range and a second 
matrix of covariances of outputs of said detectors over a second range; 

(c) second circuitry coupled to said first circuitry and connected to form an 
estimation matrix from said first and second matrices; 

(c) third circuitry coupled to said second circuitry and connected to find an 
eigenvector of said estimation matrix; and 

(c) fourth circuitry coupled to said third circuitry and connected to combine 
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said outputs relatively weighted according to the components of said 
eigenvector. 

9. The receiver of claim 8, wherein: 

(a) said first, second, third, and fourth circuitry include a programmed 
processor. 
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ABSTRACT OF THE DISCLOSURE 

Multipath relative channel estimations for weighting maximal ratio 
combining of RAKE detectors in wireless communication systems uses maximal 
eigenvectors of covariance matrices of path signals. Estimates for close-in and 
outlying sets of symbols provides linear time change channel estimation. 
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